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ABSTRACT 

We present results of the analysis of 70 RR Lyrae stars located in the bar of the Large Mag¬ 
ellanic Cloud (LMC). Combining spectroscopically determined metallicity of these stars from 
the literature with precise periods from the OGLE III catalogue and multi-epoch K s photometry 
from the VISTA survey of the Magellanic Clouds system (VMC ), we derive a new near-infrared 
period-luminosity-metallicity (PLk s Z) relation for RR Lyrae variables. In order to tit the rela¬ 
tion we use a fitting method developed specifically for this study. The zero-point of the relation 
is estimated in two different ways: by assuming the value of the distance to the LMC and by 
using Hubble Space Telescope (HST) parallaxes of five RR Lyrae stars in the Milky Way (MW). 
The difference in distance moduli derived by applying these two approaches is ~ 0.2 mag. To 
investigate this point further we derive the PLk s Z relation based on 23 MW RR Lyrae stars 
which had been analysed in Baade-Wesselink studies. We compared the derived PLk s Z rela¬ 
tions for RR Lyrae stars in the MW and LMC. Slopes and zero-points arc different, but still 
consistent within the errors. The shallow slope of the metallicity term is confirmed by both 
LMC and MW variables. 

The astrometric space mission Gaia is expected to provide a huge contribution to the de¬ 
termination of the RR Lyrae PLk s Z relation, however, calculating an absolute magnitude from 
the trigonometric parallax of each star and fitting a PLk s Z relation directly to period and ab¬ 
solute magnitude leads to biased results. We present a tool to achieve an unbiased solution by 
modelling the data and inferring the slope and zero-point of the relation via statistical methods. 

Subject headings: stars: variables: RR Lyrae - galaxies: Magellanic Clouds - (cosmology): 
distance scale - methods: data analysis - stars: statistics - astrometry 


1. Introduction 

RR Lyrae stars arc radially pulsating variables connected with low-mass helium-burning stars on the 
horizontal branch (HB) of the colour-magnitude diagram (CMD). These objects are Population II stars, 
which are abundant in globular clusters and in the halos of galaxies. RR Lyrae stars arc a perfect tool for 
studying the age, formation and structure of their parent stellar system. Moreover, they arc widely used for 
the determination of distances in the Milky Way (MW) and to Local Group galaxies. 


'Based on observations made with VISTA at ESO under programme ID 179.B-2003. 
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RR Lyrae stars are primary distance indicators because of the existence of a narrow luminosity-metallicity 
(My - [Fe/H]) relation in the visual band and of period-luminosity-metallicity (PLZ) relations in the in¬ 


frared passbands. The near-infrared PLk s Z relation of RR Lyrae stars was originally discovered by Long- 


more et al. (1986 1, and later was the subject of study by many different authors (e.g., Bono et al. 2003 


Catelan et al.|2004] Del Principe et al.|2006[ Sollima et al.|2006~ Sollima et al.|2008| |Borissova et al.[2 009 


Coppola et al. 2011 Ripepi et al. 2012a). The near-infrared PLk s Z relation has many advantages in compar¬ 


ison with the visual My - [Fe/H] relation. First of all, the luminosity of RR Lyrae stars in the K s passband 
is less dependent on metallicity and interstellar extinction (Ak. ~ O.lAy). Furthermore, light curves of 
RR Lyrae stars in the K & band have smaller amplitudes and arc more symmetrical than optical light curves, 
making the determination of the mean K s magnitudes easier and more precise. 

In order to calibrate the PLk s Z relation a large sample of RR Lyrae stars is required, spanning a wide 
range of metallicities, for which accurate mean K s and [Fe/H] measurements are available. We have selected 
70 RR Lyrae variables in the Large Magellanic Cloud (LMC) with spectroscopically determined metallicities 
in the range of -2.06 < [Fe/H] < -0.63 dex (Gratton et alT] 2004 >. All of them have counterparts in 
the OGLE III catalogue (Soszyriski et al. 20091, therefore very precise periods are available. In order to 
increase the accuracy of the determination of mean K s magnitudes, multi-epoch photometry is needed. For 
this reason we are using data from the near-infrared VISTA survey of the Magellanic Clouds System (VMC, 


Cioni et al. 2011), which is performing k\,-band observations of the whole Magellanic System in 12 (or 
more) epochs, while in many previous studies only single-epoch photometry from the Two Micron All-Sky 


Survey (2MASS, Cutri et al. 2003) was used. To fit the PL[^Z relation we apply a fitting approach developed 
for the current study. This method takes into account errors in two dimensions, the intrinsic dispersion of 
the data and the possibility of inaccuracy in the formal error estimates. 

One main issue in the determination of distances with the RR Lyrae PLk s Z relation is the calibration 
of the zero-point. Trigonometric parallaxes remain the only direct method of determining distances to 
astronomical sources, free of any assumptions (such as, for instance, the distance to the LMC, etc.) and 
hence calibrating the PLk s Z zero-point. However, reasonably well estimated parallaxes exist, so far, only 
for five RR Lyrae variables in the MW observed by Benedict et al.|( 2011| > with the Hubble Space Telescope 
Fine Guidance Sensor {HST/ FGS). In this study we use both a global estimate of the LMC distance and 
the HST parallaxes in order to calibrate the zero-point of our PLk s Z relation based on LMC RR Lyrae 
stars. Furthermore, to check whether the RR Lyrae PLk s Z relation is universal and could thus be applied 
to measure distances in the MW and to other galaxies, we analyse a sample of 23 MW RR Lyrae stars, for 
which absolute magnitudes in the K and V passbands are available from the Baade-Wesselink studies (e.g., 
Femley et al. ||1998b[ and references therein). Based on these absolute magnitudes and applying our fitting 
approach we fit the RR Lyrae PLk s Z relation. Then we compare the PLk s Z relations derived for RR Lyrae 
stars in the MW and in the LMC. 


Gaia, the European Space Agency (ESA) cornerstone mission launched in December 2013, is expected 
to provide a great contribution to the determination of the RR Lyrae PLk s Z relation and to the definition of its 
zero-point in particular. The satellite is designed to produce the most precise three-dimensional (3D) map of 
the MW to date (Perryman et al7|2001 [) by measuring parallaxes of over one billion stars during its five-year 
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mission, among which are thousands of RR Lyrae variables. In the current study we present a method which 
avoids the problems of the non-linear transformation of trigonometric parallaxes (and negative parallaxes) 
to absolute magnitudes, and apply this method to fit the PLk s Z relation of the 23 MW RR Lyrae stars, based 
on simulated Gaia parallaxes. 

In Section 2 we provide information about the 70 RR Lyrae stars in the LMC that form the basis of 
the present study. In Section 3 we present our method and results of fitting the RR Lyrae PLk s Z relation in 
the LMC and in the MW. In Section 4 we present the method to fit the PLk s Z relation with simulated Gaia 
parallaxes and apply this method to the 23 MW RR Lyrae stars analysed in Section 3. Section 5 provides 
a summary of the results. In the Appendix sections we present a detailed description of the fitting method 
which was developed for this study (Appendix A) and a compilation of metal abundances for the MW RR 
Lyrae stars (Appendix B). 


2. Data 


Optical photometry for the LMC RR Lyrae stars discussed in this paper was obtained by Clementini 
et al. ( 20031 and Di Fabrizio et al. (2O05j) using the Danish 1.54 meter telescope in La Silla, Chile. Two 
different sky positions, hereafter called fields A and B were observed. Both are located close to the bar of the 
LMC ( Clementini et al.|2003[|Di Fabrizio et al.|2005| ). As a result, accurate B, V and I light curves tied to the 
Johnson-Cousins standard system and pulsation characteristics (period, epoch of maximum light, amplitudes 
and mean magnitudes) for 125 RR Lyrae stars were obtained (Di Fabrizio et al. 20051. Low-resolution spec¬ 
tra for 98 of these variables were collected by Gratton et al. ( j2004 i using the FOcal Reducer/low dispersion 
Spectrograph (FORS1) instrument mounted at the ESO VLT. They were used to derive metal abundances 
for individual stars by comparing the strength of the Ca II K line with that of the H lines (Preston |jl959 1. 
For the calibration of the method, four clusters with metallicity in the range [-2.06; -1.26] dex were used. 


According to Gratton et al. (2004), the obtained metallicities are tied to a scale, which is, on average, 0.06 


dex more metal-rich than the Zinn & West (1984) metallicity scale. 


We cross-matched the sample of 98 RR Lyrae variables with known metallicities against the catalogue 
of RR Lyrae stars observed by the OGLE III survey ( [Soszyriski et al. 20091. The OGLE III catalogue 
contains information about the position, photometric and pulsation properties of 24906 RR Lyrae stars in 
the LMC. We found that, respectively, 94, 2 and 2 objects are cross-identified with sources in the OGLE III 
catalogue within a pairing radius of 1", 3" and 7". The 2 stars with a counterpart at more than 5" are OGLE- 
LMC-RRLYR-10345 and OGLE-LMC-RRLYR-10509; for these two objects we checked both the OGLE III 


finding charts and Gratton et al. (2004) Figure 5 (field Bl) in order to understand if they are affected by any 


problem. Star OGLE-LMC-RRLYR-10345 is an isolated slightly elongated star without any clear blending 
problem, while star OGLE-LMC-RRLYR-10509 is very close to another source possibly making more 
difficult to locate accurately the star center. Considering that Gratton et al. ; ( 2004 1 and OGLE III periods for 
these 2 stars agree within 0.5%, we kept these stars in our sample. 


We compared the periods of the 98 RR Lyrae stars provided by Di Fabrizio et al. (20051 and those in 
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the OGLE III catalogue (Soszynski et al .|2009 1. For 96 objects the periods agree within ~ 2%, while for 
two objects periods differ significantly. For star A6332 the difference is of ~ 25% and for star A5148 it is 
of ~ 37% (star identifications are from [Di Fabrizio et al.]j2005} >. Moreover, star A5148 has been classified 
as a first-overtone RR Lyrae star (RRc) in the OGLE III catalogue, and as a fundamental-mode RR Lyrae 


(RRab) by Di Fabrizio et al. (2005). Since accurately estimated periods and classifications play a key role 


in the current study, we discarded these two objects from the following analysis. 

Seven objects (B2811, B4008, B3625, B2517, A2623, A2119, A10360) from the sample are classified 


as RRc by Di Fabrizio et al. (2005) and as second-overtone RR Lyrae star (RRe) in the OGLE III catalogue. 


We removed them from our analysis because of the uncertain classification. Furthermore, since one of the 
main purposes of the current research is to study the PLk s Z relation of RR Lyrae stars of ab- and c-types we 


discarded seven objects, which were classified as double-mode RR Lyrae stars (RRd) by Di Fabrizio et al. 


(2005): A7137, A8654, A3155, A4420, B7467, B6470 and B3347. This left us with a final sample of 61 


RRab and 21 RRc stars, which all have a counterpart in the OGLE III catalogue. The period search for the 
RR Lyrae stars in the OGLE III catalogue was performed using an algorithm based on the Fourier analysis 
of the light curves (Soszynski et al. |2009 1. The uncertainties in the OGLE III periods for the 82 RR Lyrae 
stars in our sample are declared to be less than 5 x 10 ~ 6 days. Therefore we used the periods provided by 
the OGLE III catalogue in order to fit the PLk s Z relation for our sample, and did not consider errors in the 
periods since they are negligible in comparison to the other uncertainties. 

In order to derive mean K s magnitudes for the RR Lyrae stars in our sample we used data from the VMC 
survey (Cioni et al. 2011). Started in 2009, the VMC survey covers a total area of 116 deg 2 in the LMC 
with 68 contiguous tiles. The survey is obtaining YJK S photometry. The V s -band photometry is taken in 
time-series mode over 12 (or more) separate epochs and each single epoch reaches a limiting K s magnitude 
~ 19.3 mag with a S/N ~ 5 (see Figure 1 of Moretti et ak| [2014 i. On the bright side, VMC is limited by 
saturation at K s ~ 11.4 mag. The majority of RR Lyrae stars in our sample are located within the VMC 
tile LMC 5_5. Observations of the tile LMC 5_5 were performed in 15 epochs taken in the period from 
2010, October 30, to 2012, January 11. For two epochs of observation the ellipticity was too high, so these 
data were not considered in the analysis. Among the remaining 13 epochs there are 11 deep and 2 shallow 
epochs. Since shallow observations were obtained in good seeing conditions their S/N was enough to detect 
the RR Lyrae stars. In the following analysis we used all 13 available epochs to fit the light curves of the 
RR Lyrae stars. PSF photometry of the time-series data for this tile was performed on the homogenised 
epoch-tile images ( Rubele et al.||2012 ) using the IRAF Daophot packages (Stetson et al. 1990). On each 
epoch-tile image the PSF model was created using 2500 stars uniformly distributed, finally the Daophot 
ALLSTAR routine was used to perform the PSF photometry on all epoch images and time-series catalogues 
were correlated within a tolerance of one arcsec. 


We have cross-matched our sample of 82 RR Lyrae stars against the PSF photometry catalogue of the 
VMC tile LMC 5_5. VMC counterparts for 71 objects were found within a pairing radius of 1". Among 
them, 70 sources have 13 epochs in the V s -band, while for one object (B4749) we have observations only 
in 6 epochs. Six data points are not enough for a reliable fit of the light curve and, consequently, for 
the robust determination of the mean K s magnitude, hence, we discarded this source from the following 
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analysis and proceeded with the 70 RR Lyrae stars, for which 13 epochs in the Af s -band exist. We derived 
the mean magnitudes of these 70 RR Lyrae stars by Fourier fitting the light curves with the GRaphical 
Analyzer of Times Series package (GRATIS, custom software developed at the Observatory of Bologna 
by R Montegriffo, see e.g. Clementini et al. 12000 1 . To fit the light curves we discarded obvious outliers. 
Nevertheless, after the cr-clipping procedure, each source still has 11 or more data points. Examples of the 
K s light curves are shown in Figure [T] 

After deriving mean magnitudes we performed the dereddening procedure. Clementini et al. (20031 
estimated reddening values of E(B — V) — 0.116 ± 0.017 and 0.086 ± 0.017 mag in LMC field A and 


B, respectively, using the method from Sturch (1966) and the colours of the edges of the instability strip 


defined by the RR Lyrae variables. Applying the coefficients from Cardelli et al. ( f 1989' 1 of A k/A y = 0.114 
and assuming a ratio of total to selective absorption of Ry = 3.1, we estimated the extinction in the 7f s -band 


as: 


A Ks = 0.35 x E(B - V) 


( 1 ) 


Table [T] summarizes the properties of the sample of 70 RR Lyrae stars which have a counterpart in the 


VMC catalogue. First and second columns give the identification of the stars in Di Fabrizio et al. (2005) and 
in the OGLE III catalogue, respectively. The table also shows coordinates and the classification of the stars 


from the OGLE III catalogue, metallicities with errors from Gratton et al. (2004) and dereddened mean K s 
magnitudes, determined with the GRATIS package, along with their errors. 





















-7- 



17.8 


1 1 I 1 1 1 1 I 1 1 1 1 I 1 ' 

A26821 P = 0.5969220 


18 - 


80 

cd 






fi 18.2 - 

73 

^ 18.4 - 






0 


0.5 1 

Phase 


1.5 




_ i i i i | i i i 1 | 1 1 1 1 | 1 i 1 i _ 



I ' ' ' I ' ' ' ' I ' I I I I I I I I 


18 

1 B4179 P = 0.3545232 


17.8 

4 A7477 P = 0.6564084 

"go 

cd 


i\\ s|| i 

cd 




a 

18.2 

s 5 S 

a 

18 

m 9 

73 

18.4 


73 

18.2 





1 1 1 1 1 1 1 1 1 1 1 1 1 



- i i i i l i i i i 1 i i i i 1 i i i i - 


0 0.5 1 1.5 2 0 0.5 1 1.5 2 

Phase Phase 


&o 

cd 


£ 


73 




17.8 

18 

18.2 


_ 1 1 1 1 i 1 1 1 


_ 



_ 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 _ 

1 B1907 P 

= 0.5818283 



17.8 

4 A9154 P = 0.6182903 4 


# tt. 

• h 

I 

(mag) 

18 


m 

9 


73 

18.2 

: : 

i l l 1 : i 

_l_1_1_1_1_1_1_1_1_L 

, - 


18.4 

“ , , , , 1 , , i , 1 i , , i 1 .... ~ 


0 0.5 1 1.5 2 

Phase 


0 0.5 1 1.5 2 

Phase 


Fig. 1.— Examples of X^-band light curves for RR Lyrae stars in our sample. Identification numbers are 


from Di Fabrizio et al.|(|2005), periods are from the OGLE III catalogue (Soszynski et al. 2009) and are 


given in days. 














Table 1. Properties of the 70 RR Lyrae stars in the bar of the LMC analyzed in this paper 


Star 

OGLE ID 

RA 

(J2000) 

DEC 

(J2000) 

Type 

[Fe/H] 

(dex) 

0"[Fe/H] 

(dex) 

P 

(days) 

(K s ,o) 

(mag) 

0-<x,,o> 

(mag) 

A28665 

OGLE-LMC-RRLYR-12944 

5:22:06.55 

-70:27:55.6 

RRc 

-0.63 

0.24 

0.3008299 

18.450 

0.046 

A7864 

OGLE-LMC-RRLYR-13857 

5:23:39.25 

-70:31:38.1 

RRc 

-1.36 

0.22 

0.3129458 

18.550 

0.055 

B4946 

OGLE-LMC-RRLYR-10621 

5:18:11.08 

-70:59:35.6 

RRc 

-1.11 

0.25 

0.3130142 

18.394 

0.054 

A2636 

OGLE-LMC-RRLYR-13548 

5:23:09.09 

-70:39:08.1 

RRc 

-1.61 

0.29 

0.3154437 

18.562 

0.045 

A8837 

OGLE-LMC-RRLYR-13326 

5:22:45.70 

-70:30:14.3 

RRc 

-1.52 

0.22 

0.3165579 

18.660 

0.089 

A8622 

OGLE-LMC-RRLYR-13164 

5:22:28.93 

-70:30:35.9 

RRc 

-1.44 

0.28 

0.3212334 

18.426 

0.032 

A7231 

OGLE-LMC-RRLYR-13680 

5:23:22.42 

-70:32:35.4 

RRc 

-1.46 

0.26 

0.3228047 

18.236 

0.051 

A2234 

OGLE-LMC-RRLYR-13479 

5:23:01.47 

-70:39:44.4 

RRc 

-1.53 

0.18 

0.3228060 

18.292 

0.044 

A4388 

OGLE-LMC-RRLYR-12614 

5:21:31.67 

-70:36:46.3 

RRc 

-1.33 

0.27 

0.3417737 

18.411 

0.048 

A10113 

OGLE-LMC-RRLYR-14046 

5:24:00.38 

-70:28:06.1 

RRc 

-1.52 

0.25 

0.3506618 

18.197 

0.036 

B6255 

OGLE-LMC-RRLYR-10111 

5:17:17.88 

-70:57:26.4 

RRc 

-1.52 

0.16 

0.3535596 

18.305 

0.038 

B4179 

OGLE-LMC-RRLYR-10142 

5:17:19.95 

-71:01:02.1 

RRc 

-1.53 

0.27 

0.3545232 

18.150 

0.034 

A8812 

OGLE-LMC-RRLYR-13150 

5:22:26.44 

-70:30:19.1 

RRc 

-1.23 

0.24 

0.3549660 

18.281 

0.036 

A26715 

OGLE-LMC-RRLYR-12593 

5:21:29.33 

-70:29:23.4 

RRc 

-1.39 

0.18 

0.3569006 

18.308 

0.032 

A2024 

OGLE-LMC-RRLYR-13572 

5:23:11.02 

-70:40:03.3 

RRc 

-1.62 

0.26 

0.3590534 

18.246 

0.043 

B6164 

OGLE-LMC-RRLYR-10612 

5:18:10.17 

-70:57:30.7 

RRc 

-1.88 

0.22 

0.3744821 

18.039 

0.045 

A27697 

OGLE-LMC-RRLYR-13012 

5:22:14.03 

-70:28:35.0 

RRc 

-1.33 

0.25 

0.3825700 

18.030 

0.023 

A19450 

OGLE-LMC-RRLYR-13841 

5:23:37.95 

-70:34:06.7 

RRab 

-0.76 

0.13 

0.3979182 

18.481 

0.071 

B7064 

OGLE-LMC-RRLYR-10708 

5:18:18.63 

-70:55:58.7 

RRc 

-2.03 

0.20 

0.4004744 

18.029 

0.043 

B6957 

OGLE-LMC-RRLYR-10702 

5:18:18.08 

-70:56:08.7 

RRc 

-1.48 

0.18 

0.4047399 

18.027 

0.043 

B23502 

OGLE-LMC-RRLYR-10509 

5:18:00.25 

-70:54:31.0 

RRab 

-1.55 

0.14 

0.4724681 

18.243 

0.076 

A3061 

OGLE-LMC-RRLYR-13704 

5:23:25.18 

-70:38:28.9 

RRab 

-1.26 

0.12 

0.4744410 

18.415 

0.041 

B10811 

OGLE-LMC-RRLYR-10684 

5:18:16.01 

-71:04:27.0 

RRab 

-1.42 

0.20 

0.4760753 

18.197 

0.036 

B3400 

OGLE-LMC-RRLYR-10072 

5:17:14.51 

-71:02:26.6 

RRab 

-1.45 

0.24 

0.4852148 

18.346 

0.092 

A7325 

OGLE-LMC-RRLYR-13855 

5:23:39.13 

-70:32:24.8 

RRab 

-1.18 

0.26 

0.4864544 

18.223 

0.046 

B3033 

OGLE-LMC-RRLYR-10659 

5:18:14.04 

-71:03:00.5 

RRab 

-1.26 

0.21 

0.4986975 

18.130 

0.066 

B2055 

OGLE-LMC-RRLYR-10108 

5:17:17.44 

-71:04:50.2 

RRab 

-1.70 

0.23 

0.5207746 

18.254 

0.074 

A26525 

OGLE-LMC-RRLYR-12811 

5:21:52.50 

-70:29:28.7 

RRab 

-1.41 

0.22 

0.5225029 

18.168 

0.053 

A7211 

OGLE-LMC-RRLYR-13092 

5:22:21.17 

-70:32:43.9 

RRab 

-1.33 

0.19 

0.5226857 

18.193 

0.041 

A2767 

OGLE-LMC-RRLYR-13634 

5:23:17.75 

-70:38:55.9 

RRab 

-1.37 

0.08 

0.5325871 

18.054 

0.036 

B24089 

OGLE-LMC-RRLYR-10345 

5:17:43.51 

-70:54:02.7 

RRab 

-1.48 

0.16 

0.5580613 

18.094 

0.069 

A8788 

OGLE-LMC-RRLYR-13678 

5:23:22.41 

-70:30:14.6 

RRab 

-1.61 

0.21 

0.5591710 

18.197 

0.036 

A6398 

OGLE-LMC-RRLYR-13294 

5:22:40.76 

-70:33:50.2 

RRab 

-1.40 

0.30 

0.5619466 

17.957 

0.027 

A7247 

OGLE-LMC-RRLYR-13708 

5:23:25.58 

-70:32:33.4 

RRab 

-1.38 

0.21 

0.5621512 

18.045 

0.049 

A25301 

OGLE-LMC-RRLYR-12638 

5:21:34.00 

-70:30:24.5 

RRab 

-1.58 

0.27 

0.5631146 

18.268 

0.051 

A15387 

OGLE-LMC-RRLYR-12603 

5:21:30.43 

-70:37:11.3 

RRab 

-1.81 

0.12 

0.5635914 

18.049 

0.052 

B22917 

OGLE-LMC-RRLYR-10713 

5:18:19.10 

-70:54:56.1 

RRab 

-1.29 

0.16 

0.5646803 

18.179 

0.054 

A9245 

OGLE-LMC-RRLYR-13536 

5:23:07.67 

-70:29:36.5 

RRab 

-1.27 

0.18 

0.5678763 

18.013 

0.035 

A12896 

OGLE-LMC-RRLYR-13330 

5:22:46.15 

-70:38:54.9 

RRab 

-1.53 

0.10 

0.5719281 

18.143 

0.035 

A7609 

OGLE-LMC-RRLYR-13941 

5:23:48.39 

-70:32:00.3 

RRab 

-1.63 

0.11 

0.5724984 

18.023 

0.046 

B7442 

OGLE-LMC-RRLYR-10082 

5:17:15.73 

-70:55:26.8 

RRab 

-1.58 

0.11 

0.5740274 

18.096 

0.048 

A25362 

OGLE-LMC-RRLYR-13848 

5:23:38.53 

-70:30:08.5 

RRab 

-1.39 

0.15 

0.5787944 

18.033 

0.052 

B1907 

OGLE-LMC-RRLYR-10638 

5:18:12.36 

-71:04:59.5 

RRab 

-1.70 

0.26 

0.5818283 

17.915 

0.036 

A4974 

OGLE-LMC-RRLYR-13372 

5:22:51.26 

-70:35:47.7 

RRab 

-1.36 

0.10 

0.5820430 

17.992 

0.054 
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3. PLk Z relation of RR Lyrae stars 
3.1. Method 


Using the dereddened mean K s magnitudes of the 70 RR Lyrae stars derived as described in Section[2] 
spectroscopically determined metallicities from Gratton et al. I(2004 1 and accurately estimated periods from 
the OGLE III catalogue (with RRc stars "fundamentalized" by adding 0.127 to the logarithm of the period) 
we can now fit the PLk s Z relation. The fit is performed using a fitting approach developed specifically for 
this work. 


Fitting a line to data is a common exercise in science. Most common approaches use Minimum- 
Least-Squares methods, however these are often based on assumptions which do not always hold for real 
observational data. The most basic methods assume that data are drawn from a thin line with errors, which 
are Gaussian, perfectly known, and exist in one axis only. These assumptions do not hold in the present 
case, as we have an unknown but potentially significant intrinsic dispersion, non-negligible errors in two 
dimensions (K s and [Fe/H]), and the possibility of inaccuracy in the formal error estimates (e.g. in the 
determination of the precision metallicity estimates). 


We therefore follow the prescription of Hogg et al. ( 2010 ), who develop a method for fitting a line 
to data which avoids the problems highlighted above by statistical modelling of the data. They present a 
method for use in two dimensions, which has been extended to three dimensions in this paper. 


The method assumes that the data is drawn from a plane of the form 


Ls(P, [Fe/H]) = A logP + B [Fe/H] + C 


( 2 ) 


where A is the slope in the logP axis, B is the slope in the metallicity [Fe/H] axis, and C is the intercept. 
We assume a uniform Gaussian intrinsic dispersion around the luminosity axis, plus the scatter caused by 
the Gaussian observational errors. The exact mathematical definition is given in Appendix [A] The method 
utilises adaptive Markov Chain Monte Carlo (MCMC) methods (Foreman-Mackey et al .|2013[ ) to evaluate 
the posterior probability density function (PDF) of each parameter, given an input dataset, and returns the 
maximum a posteriori probability (MAP) estimate of each parameter, the formal error estimate, and the full 
posterior PDF. The formal error estimate is obtained from the 16% and 84% qualifies of the posterior PDF of 
the parameters, which give the 1 c r formal error estimate assuming that the posterior PDF is approximately 
Normal. The free fit parameters are: the slope in logP, the slope in metallicity, the zero-point, and the 
intrinsic dispersion perpendicular to the magnitude axis. 

By applying this method we found the following relation between period, metallicity and mean apparent 
K s magnitude: 


£ s ,o = (-2.73 ± 0.25)logP + (0.03 ± 0.07)[Fe/H] 
+ (17.43 ±0.01) 


( 3 ) 
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Table 1—Continued 


Star 

OGLE ID 

RA 

(J2000) 

DEC 

(J2000) 

Type 

[Fe/H] 

(dex) 

0"[Fe/HJ 

(dex) 

P 

(days) 

<*s,0> 

(mag) 

(mag) 

B6798 

OGLE-LMC-RRLYR-10044 

5:17:11.37 

-70:56:32.6 

RRab 

-1.40 

0.23 

0.5822610 

17.910 

0.114 

B14449 

OGLE-LMC-RRLYR-09999 

5:17:05.37 

-71:01:40.9 

RRab 

-1.70 

0.13 

0.5822854 

18.118 

0.071 

A9494 

OGLE-LMC-RRLYR-13354 

5:22:49.26 

-70:29:13.5 

RRab 

-1.69 

0.28 

0.5844615 

17.874 

0.035 

A18314 

OGLE-LMC-RRLYR-13353 

5:22:49.13 

-70:34:59.2 

RRab 

-1.42 

0.18 

0.5875708 

18.093 

0.030 

A10487 

OGLE-LMC-RRLYR-13126 

5:22:24.61 

-70:27:40.6 

RRab 

-1.49 

0.11 

0.5909585 

18.030 

0.016 

A10214 

OGLE-LMC-RRLYR-12609 

5:21:31.14 

-70:28:12.0 

RRab 

-1.48 

0.12 

0.5918196 

17.904 

0.065 

A28066 

OGLE-LMC-RRLYR-13765 

5:23:30.10 

-70:28:11.0 

RRab 

-1.44 

0.17 

0.5959296 

18.007 

0.060 

A26821 

OGLE-LMC-RRLYR-12831 

5:21:53.95 

-70:29:17.5 

RRab 

-1.37 

0.13 

0.5969220 

18.130 

0.047 

B2249 

OGLE-LMC-RRLYR-10061 

5:17:13.06 

-71:04:27.1 

RRab 

-1.56 

0.15 

0.6030630 

17.999 

0.050 

A16249 

OGLE-LMC-RRLYR-12960 

5:22:08.27 

-70:36:31.0 

RRab 

-1.87 

0.12 

0.6067385 

18.060 

0.045 

A4933 

OGLE-LMC-RRLYR-13175 

5:22:30.05 

-70:35:53.7 

RRab 

-1.48 

0.12 

0.6134920 

17.768 

0.027 

A7734 

OGLE-LMC-RRLYR-12956 

5:22:07.86 

-70:31:59.8 

RRab 

-1.40 

0.15 

0.6149615 

17.888 

0.036 

A2525 

OGLE-LMC-RRLYR-13788 

5:23:32.45 

-70:39:15.3 

RRab 

-2.06 

0.14 

0.6161452 

17.964 

0.051 

A9154 

OGLE-LMC-RRLYR-13494 

5:23:02.93 

-70:29:44.6 

RRab 

-1.66 

0.14 

0.6182903 

17.972 

0.029 

B1408 

OGLE-LMC-RRLYR-10067 

5:17:13.84 

-71:06:06.9 

RRab 

-1.70 

0.11 

0.6297088 

18.021 

0.013 

A5589 

OGLE-LMC-RRLYR-12968 

5:22:09.60 

-70:35:02.5 

RRab 

-1.60 

0.13 

0.6375745 

17.948 

0.035 

A7468 

OGLE-LMC-RRLYR-13176 

5:22:30.06 

-70:32:20.6 

RRab 

-1.55 

0.11 

0.6386908 

18.043 

0.041 

A25510 

OGLE-LMC-RRLYR-13002 

5:22:13.43 

-70:30:11.4 

RRab 

-1.72 

0.11 

0.6495506 

17.713 

0.038 

A8720 

OGLE-LMC-RRLYR-13956 

5:23:50.19 

-70:30:16.7 

RRab 

-1.88 

0.34 

0.6508174 

17.847 

0.037 

B7063 

OGLE-LMC-RRLYR-10973 

5:18:44.05 

-70:55:55.8 

RRab 

-1.49 

0.14 

0.6548698 

17.867 

0.023 

B7620 

OGLE-LMC-RRLYR-10541 

5:18:03.58 

-70:55:03.1 

RRab 

-2.05 

0.12 

0.6561602 

17.689 

0.034 

A7477 

OGLE-LMC-RRLYR-14068 

5:24:02.97 

-70:32:08.6 

RRab 

-1.67 

0.28 

0.6564084 

17.955 

0.017 

A28293 

OGLE-LMC-RRLYR-12758 

5:21:46.13 

-70:28:13.3 

RRab 

-1.74 

0.10 

0.6602890 

17.979 

0.050 

A6426 

OGLE-LMC-RRLYR-13196 

5:22:32.51 

-70:33:48.7 

RRab 

-1.59 

0.09 

0.6622400 

17.868 

0.038 

A3948 

OGLE-LMC-RRLYR-13285 

5:22:40.40 

-70:37:17.0 

RRab 

-1.46 

0.12 

0.6623845 

17.944 

0.036 

A8094 

OGLE-LMC-RRLYR-13306 

5:22:43.06 

-70:31:23.7 

RRab 

-1.83 

0.12 

0.7420663 

17.890 

0.033 


Note. — Columns report: 1) Star identification from |Di Fabrizio et al.|j2005) : 2) Identification from the OGLE III catalogue (Soszynski et 
al.|2009|: 3) Right ascension (OGLE); 4) Declination (OGLE); 5) Type; 6) Metallicity from|Gratton et al.||2004|; 7) Metallicity error |( I ml t on 
et al.|2004) ; 8) Period (OGLE); 9) Dereddened mean K s magnitude from the VMC data, determined from the analysis of the light curve with 
GRATIS; 10) Error of the mean K$ magnitude. 
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The intrinsic dispersion of the relation is found to be 0.01 mag. The RMS deviation of the data around 
the relation, neglecting the intrinsic dispersion, is 0.1 mag. Since the reddening in the K s band is negligible 
we suggest that the effects of the LMC depth cause the intrinsic dispersion of the relation. The left panel of 
Figure [2] presents the PLk s Z relation (Equation[3]) of the 70 LMC RR Lyrae stars in the period-luminosity- 
metallicity space, whereas the right panels show the projection of the PLk s Z on the log(P) - K s (top-right 
panel) and K s - [Fe/H] (bottom-right panel) planes. The grey lines in the figure are lines of equal metallicity 
(top-right) or equal period (bottom-right). The method finds the relation (values of A, B, and C for the 
relation K s = AlogP+B [Fe/H]+C) in the three dimensions (logP, K s , and [Fe/H]). Each of the grey lines in 
the top-right plot are therefore K s =AlogP+B [Fe/H]+C for the full range of periods, at the metallicity of 
each star (one line per star). Thus, by following the line up and down it is seen how K s changes with period 
at some specific metallicity. The lines do not always cross the points on the diagram because the line is the 
result of the fit, and the points are affected by errors and intrinsic dispersion so may be above or below the 
fit. In the bottom-right plot the lines are K s = AlogP+B [Fe/H]+C for the full range of metallicity with logP 
taken from each star. 


It is worth noting that we find a very small dependence of the K s magnitude on metallicity. However, 
the metal abundance range spanned by the adopted sample does not reach the highest values (up to solar and 
supersolar) observed in the MW bulge and disk RR Lyrae populations. In order to study the effect of the 
adopted range of metallicities on the slope of the PLk s Z relation we derived this relation also for MW RR 


Lyrae stars. We discuss the results in Section 3.3 


3.2. Zero-point of the PLk s Z relation in the LMC 

To use the derived PLk s Z relation for determining distances it is necessary to calibrate its zero-point. 
This can be done in a number of different ways. In this paper we follow two different approaches: the 
first one is based on adopting a value for the distance of the LMC; in the second one we use the absolute 
magnitudes of Galactic RR Lyrae stars for which trigonometric parallaxes have been measured with the 
HST/ FGS. Both approaches have their advantages and disadvantages, we discuss them in the following 
sections. 


3.2.1. Zero-point based on the LMC distance 


The LMC is widely considered the first rung of the cosmic distance ladder as it contains a large number 
of different distance indicators, such as Cepheids, RR Lyrae variables, eclipsing binaries (EBs), red giant 
branch (RGB) stars, etc., allowing the galaxy distance to be determined with several independent techniques. 
Figure 8 of Benedict et al. (2002) shows an impressive summary of LMC distance moduli published during 
the ten years from 1992 to 2001. Values from 18.1 to 18.8 mag were reported in the literature, with those 
smaller than 18.5 mag supporting the so-called "short" scale, and those larger than 18.5 mag, the "long" 
one. In more recent years the dramatic progress in the calibration of the different distance indicators has 
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led the dispersion in LMC distance moduli to shrink significantly. Extreme values such as those listed in 
Benedict et ah] (2002 1 are not very often seen in the recent literature (Clementini 20081. Still a general 
consensus on the LMC distance has not been fully reached yet. Moreover, there have been significant 
concerns about a possible "publication bias" affecting the distance to the LMC ( Schaeferj2008[ Rubele et al. 
2012| Wafker||2012| ). In particular, Schaefer (20081 claimed that from 2002 to 2007 June, 31 independent 
papers reported new measurements of the distance of the LMC, and the new values clustered tightly around 
the value (m - M )o = 18.5 ± 0.1 mag, adopted by the HST Key Project on the extragalactic distance scale 
( Lreedman et ah|200T i. Schaefer (2008) considered the effects of the "publication bias" to be the most likely 
cause of the clustering of LMC distance measurements. 

A number of studies on the compilation of distances to the LMC as derived from different distance indi¬ 
cators can be found in the literature of the last 15 years (e.g., Gibson j2000 ; Benedict et ah ;2()02| Clementini 
et al.|2003 ; Schaefer|2008t de Grijs et al.|2014 ). Clementini et ah (2003 1 analysed the distance to the LMC 
measured using Population I and Population II standard candles and showed that all distance determina¬ 
tions converge within lcr error on a distance modulus (m - M )o = 18.515 ± 0.085 mag. The most recent 


compilation of LMC distance moduli is that of de Grijs et ah (2014) who compiled 233 separate distance 
determinations, published from 1990 March until 2013 December, and concluded that the canonical distance 


modulus of (m - M )o = 18.49 ± 0.09 mag may be used for all practical purposes. The compilation of de Grijs 


et ah (2014) includes the distance modulus of (m - M )o = 18.46 ± 0.03 estimated by Ripepi et ah (2012b) 
using LMC classical Cepheids observed by the VMC survey, and the recent determination of direct distances 


to eight long-period EBs in the LMC by Pietrzyhski et ah (2013), which is claimed to be accurate to within 
~ 2%: Dlmc = 49.97 ±0.19 (stat) ± 1.11 (syst) kpc, corresponding to the distance modulus (in - M) o= 
18.493 ± 0.008 (stat) ± 0.047 (syst) mag. Lurthermore, the model fitting of the light curves of different 
classes of pulsating stars in the LMC, also based on different samples and hydrodynamic al codes, provides 
values consistent with 18.5 mag (see Bono et al.|2002 Marconi & Clementini|2005| |Keller & Wood|2002| 
2006 McNamara et al.|2007| ). 


The RR Lyrae stars in our sample are located in a relatively small area close to the center of the LMC 
bar. Neglecting depth/projection effects they can be considered as being all at the same distance from us 
and close to late-type EBs, which are all located relatively close to the barycenter of the LMC as derived 
by Pietrzyhski et ah (j20131. Therefore, in the following we adopt for the distance modulus of the LMC 


the value published by Pietrzyhski et ah (20131. We subtracted this value from the dereddened mean K s 
apparent magnitudes of our 70 RR Lyrae stars and derived absolute magnitudes in the K s band ( Mk s ). 


Then by applying the technique described in Section 3.1 we derived the relation between K s -band absolute 
magnitudes, periods and metallicities, with the zero-point entirely based on the distance to the LMC by 
Pietrzyhski et ah (20131: 


M Ks = (-2.73 ± 0.25)logP + (0.03 ± 0.07)[Fe/H] 
- (1.06 ±0.01) 


( 4 ) 
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3.2.2. Zero-point based on trigonometric parallaxes of Galactic RR Lyrae stars 


In order to obtain an estimate of the PLk s Z relation zero-point which is independent of the distance to 
the LMC and, in turn, be able to measure the distance to this galaxy from the RR Lyrae PLk s Z relation, it is 
necessary to know the absolute magnitude of the RR Lyrae stars with good accuracy. Trigonometric paral¬ 


laxes remain the only direct method to measure distances and hence derive absolute magnitudes. Benedict et 


al. (2011) derived absolute trigonometric parallaxes for five Galactic RR Lyrae stars (RZ Cep, XZ Cyg, SU 


Dra, RR Lyr and UV Oct) with the HST/ FGS. With these parallaxes the authors estimated absolute magni¬ 


tudes in the K s and V passbands, corrected for interstellar extinction and Lutz-Kelker-Hanson bias (Lutz & 


Kelker|1973 Hanson|1979| ). Absolute magnitudes in the X s -band, periods and metallicities from Benedict et 
al. (20111, and the slopes of the relation derived in Equation [3] were used in order to determine a zero-point 


from each of these five MW RR Lyrae stars. The metallicities in Benedict et al. (2011) are in the Zinn & 
West metallicity scale and were converted to the metallicity scale in Gratton et al. (12004 ) by adding 0.06 
dex (see Section [2]). The logarithm of the period of the RRc star RZ Cep was "fundamentalized" by adding 
0.127. Then we calculated the weighted mean of the five zero-points, this corresponds to: -1.27 ±0.08 mag. 
The relation between absolute magnitudes, periods and metallicities with the zero-point based on the five 
MW RR Lyrae stars from Benedict et al. (2011[) is: 


M Ks = (-2.73 ± 0.25)logP + (0.03 ± 0.07)[Fe/H] 
- (1.27 ±0.08) 


(5) 


A recent analysis (Monson 2015, private communication) shows that there is likely a typo in Benedict 


et al. (20111 parallax for the RR Lyrae star RZ Cep. Hence, we excluded this star from the sample and 


derived the zero-point based on parallaxes of remaining four RR Lyrae stars (XZ Cyg, UV Oct, SU Dra and 
RR Lyr): 


M Ks = (-2.73 ± 0.25)logP + (0.03 ± 0.07)[Fe/H] 

- (1.25 ±0.06) (6) 

Situation improves, however there is still a difference of ~ 0.2 mag between the two zero-point obtained 
based on the distance to the LMC (Eq. [4]l and the one based on the HST parallaxes of four RR Fyrae stars 
(Eq.[6]). In fact, if we apply our PFk s Z relation with zero-point calibrated on Benedict et al.]( j2011j ) parallaxes 
(Eq. [6]) to determine the absolute magnitudes of the 70 RR Lyrae stars in our sample, we obtain a distance 
modulus for the LMC of (m - M)q - 18.68 ± 0.10 mag. This distance modulus is about 0.2 mag longer than 
the widely adopted value of (m - M) o = 18.5 mag. 


There arc a number of possible explanations for the discrepancy between zero-points. First of all, we 
should remember that Pietrzynski et al. (j2013)> results have been called into question by Schaefer (20131 






















































-14- 


who, in addition to concerns regarding possible bandwagon effects, also pointed out that Pietrzynski et al. 


(2013) distance to the LMC differs significantly from the average distance inferred from four hot, early-type 
EBs, D=47.1±1.4 kpc ((m - M)q - 18.365 ± 0.065 mag), published by Guinan et al. (1998), [Fitzpatrick 


et al.j (2002 20031 , and Ribas et alj (2002 1 . Furthermore, in using the late-type EBs to calibrate the RR 
Lyrae PLk s Z relation we have implicitly assumed that RR Lyrae stars and EBs are at same distance from us. 
ffowever, when pushing for distance comparisons at a few percent level the effects of sample size, spatial 
distribution, depth and geometric projection become important and properly accounting for the internal 
structure of the LMC may become necessary. The RR Lyrae stars in our sample could be distributed along 


the whole depth of the LMC. Moreover, RR Lyrae stars and EBs from Pietrzynski et al. (2013) could reside in 
different sub-structures of the LMC, which could be the reason for the systematic error in the determination 
of the zero-point (see e.g. Moretti et al. [2014| for different features of the LMC structure traced by classical 
Cepheids, RR Lyrae stars and hot EBs). 

On the other hand, when calibrating the zero-point by applying parallaxes of the four MW RR Lyrae 
stars by Benedict et al. (20111 we implicitly assumed that the PLk s Z relation is the same in the MW and 


in the LMC, which may not be true (see Subsection 3.31. We may also wonder whether there might be 
unknown systematic errors affecting Benedict et al.’s parallaxes. These come from HST fields, which provide 
relative and not absolute trigonometric parallaxes. Absolute parallaxes of the reference stars in each field are 
estimated via a complex procedure of fitting the spectral type and luminosity class of each star. A general 
formal error of 0.5 mas is applied to the absolute parallax of the reference stars, equal for all stars in all 
fields, and without justification. This could result in miscalculated estimates of the precision of the final 
absolute parallax measurements of the four RR Lyrae stars. The Lutz-Kelker bias is corrected a posteriori. 


In this respect it is worth of notice that, according to van Leeuwen (2007), Hipparcos parallax of RR Lyrae 
itself, the only RR Lyrae variable for which the satellite measured the parallax with high precision (± 0.64 
mas), is about 0.31 mas smaller than Benedict et al. ( 201 l| )’s parallax for the star, although consistent with it 
within the errors, hence, the corresponding distance modulus is about 0.17 mag longer. In any case, a great 
contribution to the determination of the zero-point of the RR Lyrae PLk s Z relation is expected from the ESA 
astrometric satellite Gaia. We discuss this topic in Section [4] 


3.3. PLk s Z relation of RR Lyrae stars in the MW 


In spite of many studies in the literature, it remains still unsettled whether the RR Lyrae PLk s Z is 
a universal relation. To investigate this matter we have derived the PLk s Z relation for RR Lyrae stars in 
the MW and compared it with the relation we have obtained in Section |3.2| for the LMC variables. To 
this end we selected 23 MW RR Lyrae stars which have their absolute magnitudes known from Baade- 
Wesselink (hereinafter B-W) studies based on near-infrared data (Jones et al. 1988; 1992[ Femley et al. 
1990 Liu & Janes|1990 ; Cacciari et al.|1992| Skillen et al.|1993[|Femley et al.|1994 and references therein) 
and metallicities determined from abundance analysis of high resolution spectra. Information about these 
23 RR Lyrae stars is presented in Table [2] Star’s coordinates in the table are from the SIMBAD database; 


periods, apparent V and K magnitudes and reddening E(B-V) are from [Femley et al. (1998a). The sample 
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contains two first-overtone mode RR Lyrae stars (namely, T Sex and TV Boo). As done for the LMC RRc 
stars, their periods were fundamentalized by adding +0.127 to the logarithm of the period. Absolute My 
and Mk magnitudes in Columns 10 and 12 were taken from the compilations of B-W results in Table 11 of 
Cacciari et ah ( jl992| and from Table 16 of Skillen et al.l ( fT993) ) for the variable stars: WY Ant, W Crt and RV 


Oct. According to Cacciari et al. (1992) the K magnitudes of the stars analyzed with the B-W method are in 


the Johnson photometric system. Following the discussion in Cacciari et al. (1992j) and [Skillen et al. (19931 
we retained only 23 of the original lists of 29 stars field RR Lyrae stars analysed with the B-W method, as 
stars which are likely to be evolved (DX Del, SU Dra, SS Leo, BB Pup and W Tuc) were discarded. We 
also discarded DPI Peg as there is suspect the star is a dwarf Cepheid (see Feast et a l.|2008| and discussion 


therein). Furthermore, following Fernley et al. (1994), original My and Mk values were revised (i) assuming 
for the p factor used to convert the observed pulsation velocity to true pulsation velocity in B-W analyses 
the value of p= 1.38, and (ii) multiple determinations of individual stars were averaged. 

Metal abundances with related errors are needed to apply our fitting approach. Several different spec¬ 
troscopic studies have targeted the stars in Table [2] In Appendix B we provide a summary of their results. 
The largest and most homogeneous samples are those by Clementini et al. (1995) and Lambert et al.|( jl996| . 
These Authors measured [Fe/H] abundances from high resolution spectra for several of the stars in Table [2] 
and provided recalibrations of the A5 index (Preston j[l959 l, from which metal abundances can be derived 
for the stars which lack abundance analysis. For sake of homogeneity and ease of use in this paper we 
adopt metallicities and metallicities errors for the MW RR Lyrae stars as they are listed, ready for use, in 


Table 21 of Clementini et al. (1995). These [Fe/H] values are the average of the Fel and Fell abundances, 


adopting log e(FeI)= 7.56 and log e(FeII)=7.50 for the sun. They are reported along with the related errors 
in columns 5 and 6 of Table [2] Metallicities in Lambert et al. ( |1996[ ) were obtained from the Fell abundance 
and adopting log e(FeI)= 7.51 for the sun. Corresponding [Fe/H] values and errors, are reported in columns 


2 and 3 of Table 4 in Appendix B. In Appendix B we also present the PLZ relation obtained using Lambert 
et al.|(T996])’s metallicities and our approach. 


By applying our fitting approach (Section |3Tj ) to the 23 MW RR Lyrae stars we derived the following 
PLkZ relation: 


M k = (-2.53 ± 0.36)logP + (0.07 ± 0.04)[Fe/H] 

- (0.95 + 0.14) (7) 

The intrinsic dispersion of the relation is found to be 0.007 mag. The RMS deviation of the data around 
the relation, neglecting the intrinsic dispersion, is 0.086 mag. It should be noted that the metallicities listed 
in Table [2]may differ slightly from the metal abundances used in the B-W analysis of these stars. However, 
this is not of great concern as the B-W based on near-infrared data is mildly affected by small changes 
in metallicity and reddening. We also point out that the rather large error of the logP term in Equation [7] is 
largely driven by the large errors (0.15-0.25 mag) in the X-band absolute magnitudes from the B-W analyses. 
This is confirmed by the exercise with Gaia simulated parallaxes we present in Section 4. 
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The slope in [Fe/H] in Equation [7] is higher than the slope obtained for the LMC RR Lyrae stars 
(Equations [4j [6]), although they are still consistent within the respective errors. Equation [7] was derived over 
a wide range of metallicities [-2.5; 0.17] dex, nevertheless the slope of the metallicity term remains rather 
small. Thus, the relatively small metallicity range spanned by the LMC variables could be not responsible 
for the negligible dependence on metallicity of the RR Lyrae PLk s Z relation in the LMC. 


The distribution of the 23 MW RR Lyrae stars in the period-luminosity-metallicity space and the pro¬ 
jections of the PLkZ relation (Equation [7]) on the log(P) -Mr and M% - [Fe/H] planes are shown in Figure[3] 
Grey lines are the same as in Figure [2] and are described in Section 3.1 


Some concern may arise since the K magnitudes of the 23 MW RR Lyrae stars used to derive Equation[7] 
are in the Johnson photometric system (see, Cacciari et a l.| 1992 ), whereas for the LMC RR Lyrae variables 
we have K s photometry in the VISTA system^ To address this issue we have reported in column 10 of 
Table [2]the average K s magnitudes in the 2MASS system of the 23 MW RR Lyrae stars as provided by Feast 


|et al.| ( |2008] ). The difference with the Johnson K average magnitudes listed in column 9 is very small (of the 
order of about 0.03 mag, on average) and definitely much smaller than individual errors in the B-W K -band 
absolute magnitudes of the MW variables (0.15-0.25 mag), or errors in the K s average apparent magnitudes 
of the LMC RR Lyrae stars (see column 10 of Table 1). Hence, we are confident that the difference in 
photometric system does not affect significantly our comparison. 


3.4. Comparison with the literature 


The near-infrared PLkZ relation of the RR Lyrae stars has been studied by several authors both from 
a theoretical and an observational point of view. Longmore et al. ( 1986] ) pioneering work was followed 
by |Liu & Janes (|1990l, Skillen et al.| (19931 and Jones et al. (1996 >. A comprehensive analysis of the IR 


properties of RR Lyrae stars was performed by Nemec et al. (1994). 


Some of the RR Lyrae PLkZ relations available in the literature are presented in Table [3] Bono et al. 
(20031 derived the semi-theoretical relation presented in the first row of Table[3] This theoretical relation has 
been derived from an extended set of RR Lyrae nonlinear hydrodynamic al models spanning a wide range 
of chemical compositions (Z from 0.0001 to 0.02, which approximately corresponds to [Fe/H] from -2.45 


to -0.15 dex). Catelan et al. (2004) presented a theoretical calibration of the RR Lyrae PLkZ relation based 


on synthetic HB models computed for several different metallicities, fully taking into account evolutionary 
effects besides the effect of chemical composition. They derived the relation: 


M k = -2.353 log P + 0.175 log Z- 0.597 


(8) 


By using Eqs. 9 and 10 in Catelan et al. ( 2004j ) and assuming [o'/Fe]~ 0.3 (e.g., Carney 1996) we transformed 
Eq.[8]into the form, presented in the second row of Table [3] 


2 The VISTA system is tied to the 2MASS photometry (Skrutskie et al. 
depending on the J - K s colour, and being of the order of 3-4 mmag for the typical J 


2006), with the difference in K s magnitude only mildly 
K s colour of RR Lyrae stars. 
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Fig. 2.— Left panel: PLk s Z relation of the 70 LMC RR Lyrae stars (Equation [3]) analyzed in the paper, in 
the period-luminosity-metallicity space. Right panels: Projections of the PLk s Z relation (Equation[3]) on the 
log(P) - K & (top) and K s - [Fe/H] (bottom) planes. Grey lines represent lines of equal metallicity (top) and 
equal period (bottom). See text for the details. Uncertainties in the K s magnitude and [Fe/H] are omitted to 
simplify the figure, but they are provided in Table [I] 


logP 






Fig. 3.— Left panel: PMkZ relation of the 23 MW RR Fyrae stars (Equation [7]) in the period-luminosity- 
metallicity space. Right panel: Projections of the PM«Z relation (Equation[7]) on the \og(P)~ M/< (top panel) 
and Mk - [Fe/H] (bottom panel) planes. Grey lines represent lines of equal metallicities (top panel) and 
periods (bottom panel). See text for the details. Uncertainties in the Mk magnitude and [Fe/H] are omitted 
to simplify the figure, but they are provided in Table [2] 
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Table 2. Properties of 23 bright RR Lyrae stars in the Milky Way 


Star 

RA 

(deg) 

DEC 

(deg) 

P 

(day) 

[Fe/H] 

(dex) 

O' [Fe/H] 

(dex) 

E(B-V) 

(mag) 

V 

(mag) 

K 

(mag) 

K s 

(mag) 

M v 

(mag) 

M k 

(mag) 

O' M v ,Mk 
(mag) 

UUCet 

1.02135 

-16.99764 

0.606081 

-1.38 

0.08 

0.015 

12.08 

10.85 

10.837 

0.62 

-0.55 

0.15 

SW And 

5.92954 

29.40101 

0.442279 

-0.06 

0.08 

0.045 

9.71 

8.54 

8.505 

0.94 

-0.07 

0.15 

RR Cet 

23.03405 

1.34173 

0.553025 

-1.38 

0.08 

0.015 

9.73 

8.56 

8.520 

0.68 

-0.42 

0.15 

X Ari 

47.12869 

10.44590 

0.651139 

-2.50 

0.08 

0.18 

9.57 

7.95 

7.941 

0.57 

-0.575 

0.15 

AR Per 

64.32165 

47.40018 

0.425549 

-0.34 

0.16 

0.31 

10.51 

8.66 

8.642 

0.87 

-0.06 

0.25 

RXEri 

72.43455 

-15.74118 

0.587246 

-1.63 

0.16 

0.03 

9.69 

8.42 

8.429 

0.66 

-0.51 

0.25 

RR Gem 

110.38971 

30.88318 

0.397316 

-0.32 

0.16 

0.075 

11.38 

10.26 

10.275 

0.89 

-0.02 

0.25 

TTLyn 

135.78245 

44.58559 

0.597438 

-1.64 

0.16 

0.015 

9.86 

8.65 

8.611 

0.65 

-0.55 

0.15 

TSex 

148.36833 

2.05732 

0.324698 

-1.20 

0.16 

0.015 

10.04 

9.18 

9.200 

0.66 

-0.07 

0.25 

RR Feo 

151.93108 

23.99176 

0.452387 

-1.37 

0.16 

0.015 

10.73 

9.70 

9.730 

0.76 

-0.16 

0.15 

WY Ant 

154.02061 

-29.72845 

0.574312 

-1.32 

0.16 

0.06 

10.87 

9.64 

9.674 

0.55 

-0.55 

0.15 

WCrt 

171.62351 

-17.91435 

0.412013 

-0.89 

0.16 

0.03 

11.54 

10.56 

10.539 

0.96 

+0.08 

0.15 

TUUMa 

172.45205 

30.06733 

0.557659 

-1.38 

0.16 

0.015 

9.82 

8.67 

8.660 

0.70 

-0.41 

0.15 

UU Vir 

182.14613 

-0.45676 

0.475606 

-0.64 

0.16 

0.015 

10.56 

9.51 

9.414 

0.80 

-0.195 

0.15 

SW Dra 

184.44429 

69.51062 

0.569670 

-0.91 

0.16 

0.015 

10.48 

9.33 

9.319 

0.68 

-0.38 

0.15 

RV Oct 

206.63230 

-84.40177 

0.571130 

-1.92 

0.16 

0.09 

10.98 

9.51 

9.526 

0.68 

-0.40 

0.15 

TV Boo 

214.15242 

42.35992 

0.312559 

-2.31 

0.16 

0.015 

10.97 

10.22 

10.248 

0.58 

-0.20 

0.25 

RS Boo 

218.38839 

31.75462 

0.377337 

-0.37 

0.16 

0.015 

10.37 

9.45 

9.507 

0.85 

+0.00 

0.25 

VY Ser 

232.75803 

1.68382 

0.714094 

-1.71 

0.08 

0.03 

10.13 

8.78 

8.826 

0.61 

-0.665 

0.25 

V445 Oph 

246.17171 

-6.54165 

0.397023 

+0.17 

0.08 

0.195 

11.05 

9.24 

9.262 

1.09 

+0.30 

0.25 

TW Her 

268.63000 

30.41048 

0.399601 

-0.58 

0.16 

0.06 

11.28 

10.22 

10.239 

0.80 

-0.07 

0.15 

AV Peg 

328.01164 

22.57483 

0.390380 

-0.03 

0.16 

0.06 

10.50 

9.36 

9.346 

1.10 

+0.14 

0.15 

RVPhe 

352.13106 

-47.45362 

0.596413 

-1.69 

0.16 

0.015 

11.94 

10.72 

10.768 

0.86 

-0.29 

0.25 


Note. — The columns report: 1) Name of the star; 2) Right Ascension (J2000) from SIMBAD database; 3) Declination (J2000) from SIMBAD 
database; 4) Period from Femley et al. (1998a); 5) Metallicity from Cleme ntini et al.| fl995 ); 6) Error in metallicity from Clemen tini et al.| <p~995 1; 
7) Reddening from Fernley et al. ( 1 998a|i; 8) V magnitude from Femley et al. (1998aI; 9) K magnitude in the Johnson system from |Femley et al.] 
(1998a) ; 10) K.. magnitude in the 2MASS system from |Feast et al,|pb08^ ; 11) Absolute magnitude in the V passband from |Femley et al.| ( |1994|) ; 
12) Absolute magnitude in the K passband obtained from B-W analyses and corrected to p=1.38 (see text for details); 13) Errors in the absolute 
V, K magnitudes. 
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Dall'Ora et al. (2004) derived an empirical relation between apparent K magnitude and period for 21 


RRab and 9 RRc stars in the LMC globular cluster Reticulum. Del Principe et al. (2006 1 obtained the 
relation between apparent K s magnitude, metallicity and period from the analysis of RR Lyrae stars in the 
Galactic globular cluster a> Cen. 


Sollima et al.| (2006l derived a PLkZ relation by analysing 538 RR Lyrae stars in 15 Galactic clusters 
and in the LMC globular cluster Reticulum. This relation spans the metallicity range -2.15 < [Fe/H] < -0.9 


dex. Mean K magnitudes were estimated by combining Two-Micron-All-Sky-Survey (2MASS, Cutri et al. 


2003) photometry and literature data. The zero-point was calibrated on RR Lyrae itself, whose distance 


modulus was derived using the HST trigonometric parallax measured for this star by Benedict et al. (2002). 


Sollima et al. (2008) presented JKH time-series photometry of RR Lyrae and derived a new zero-point of 


Sollima et al. ( j2006j )’s PLkZ relation. 


Borissova et al. (2009) presented near-infrared K s photometry and spectroscopically measured metal¬ 


licity for a sample of 50 field RR Lyrae stars in inner regions of the LMC. These authors had 5 measurements 


in the K s passband for most of the stars in their sample and used templates from Jones et al. (1996) to fit the 
light curves and derive the mean K s magnitudes. To improve statistics they added to their sample LMC RR 
Lyrae stars from Szewczyk et al. ( 2008) dataset, and derived the PLk s Z relation based on the total sample 


of 107 LMC variables. The zero-point was calculated using Sollima et al. (2008 )’s mean K magnitude, the 
reddening and the trigonometric parallax of RR Lyrae. 


Benedict et al. ( |2011[ ) recalibrated ah the literature relations listed in Table[3] but Catelan et al. ( 2004| ’s 
one, by fitting to equations in the form: a(logP + 0.28) + Z?([Fe/H ] + 1.58) + ZP, the Lutz-Kelker-Hanson- 
corrected absolute magnitudes of the five MW RR Lyrae stars for which HST parallaxes are available. Since 
there are concerns about Benedict et al. (201T| ) RR Lyrae star (namely RZ Cep) we transformed the PLkZ 
relations of the LMC and MW RR Lyrae stars derived in this paper (Equations [4j [6]) to the form adopted 
by Benedict et al. ( 2011 1 and determined their zero-points on the basis of Benedict et al. ( 20Tl) ’s HST 
parallaxes but excluding RZ Cep. The zero-points based on Benedict et al. ( 2011 1 ) parallaxes are presented 
in column 5 of Table [3j whereas the zero-point of our LMC PLk s Z relation calculated by assuming the 


distance to the LMC (Subsection 3.2) and the zero-point of our MW PLkZ relation based on the B-W 
studies are (Section [33j ) are presented in Column 4. 

The slope in period of the RR Lyrae PLkZ relations (Column 2 of Table [3]) differs significantly in 
the various studies. The slope we derived for the LMC RR Lyrae stars is in excellent agreement with that 


derived by Del Principe et al. (2006), whereas the slope of the MW RR Lyrae PLkZ is in good agreement 


with that derived by Solhma et al. (2006 20081. 


The dependence on metallicity of the PLkZ relations (Column 3 of Table [3]) also varies among different 
studies and generally is larger in the theoretical and semi-theoretical relations. The comparison of the metal¬ 
licity dependence in the different empirical relations is complicated by the inhomogeneity of the metallicity 
scales adopted in these studies. Metallicities in Del Principe et al.J (2006) are in the |Zinn & West (1984 > 


scale, whereas in Sollima et al. (2006 2008) are in the Carretta & Grafton (19971 scale. In the current study 


for the LMC RR Lyrae stars we used the metallicity scale defined by Gratton et al.:( 2004j ) which is also the 
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scale adopted by Borissova et al.:( 2009j ). As discussed in Gratton et al. (2004) this scale is systematically 
0.06 dex higher than Zinn & West scale. This difference is small and systematic, hence should not affect 
the results of this comparison. Finally, for the MW RR Lyrae stars we used the metallicities measured by 


Clementini et al. (1995). Because the spectroscopic [Fe/H] values in Clementini et al. (1995) are derived 


from high dispersion spectra analyzed using standard reduction procedures, the derived metallicities are on 
the scale of the high dispersion spectra (i.e., the Carretta et a l.|2009 scale) and could be transformed to Zinn 
& West scale using the relations provided in Carretta et al. ( 2009j ). 

The slope in metallicity of the PLkZ relation based on the LMC RR Lyrae stars is the smallest among 
the various studies listed in Tableland it is close to Borissova et alj(2009l’s slope. This is consistent with 
the two studies both involving LMC variables and using exactly the same metallicity scale. The slope in 
metallicity we found for the MW RR Lyrae stars is larger than that of the LMC RR Lyrae stars and, in spite 
of the difference in metallicity scales, it is very close to the slope obtained by Sollima et al. (2006) for RR 
Lyrae stars in globular clusters, ffowever, taken at face value, the metallicity slopes of the empirical relations 
in Table [3] appear to be all rather small and in agreement to each other within the relative uncertainties, thus 
generally suggesting a mild dependence the RR Lyrae PLkZ independently on the specific environment. 


4. Gaia observation of RR Lyrae stars in the Milky Way 


The Gaia astrometric satellite will revolutionise many fields of astronomy (Perryman et a l.|2001 I. Of 
particular importance will be its catalogue of trigonometric parallaxes for more than one billion stars, with 
astrometric precision down to pas level. Due to Gaia’s constant observation of the sky over the five year 
nominal mission, Gaia will repeatedly observe all stars brighter than its limiting magnitude, with an average 
of 70 observations per star. This will also make it possible for Gaia to discover and characterise many types 
of variables, including RR Lyrae stars and Cepheids. 


Gaia is observing in the broad visual band G (Jordi et al. !2() 10 > for its astrometric measurements, and is 
therefore not ideal for characterising the RR Lyrae PLZ relation, which exists only in the infrared passbands. 
ffowever, since Gaia will provide accurate parallaxes for an expected tens of thousands of MW RR Lyrae 
stars, it could serve as a perfect tool for the determination of the zero-point of the PLk s Z relation through 
a combination with external datasets. As it was discussed in Subsection 3.2.2[ the current largest limiting 
factor in zero-point calibration of PLk s Z and My - [Fe/H] relations is the lack of a reliable and statistically 
significant sample of parallax measurements. The current state of the art is the sample of five RR Lyrae 


parallaxes from Benedict et al. (2011) using the HST. Gaia will improve this situation by several orders of 
magnitude in both precision and numbers of objects. Moreover, the distance to the LMC will be determined 
through the combination of Gaia parallaxes for a large sample of LMC bright stars, hence, a zero-point of 
the PLk s Z relation based on the distance to the LMC will be obtained with a high precision. 
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Table 3. PLk Z relations from the literature 


Relation 


a 


b 


ZP 1 from the original relation 


ZP 2 from 


Benedict et al. 


(2011 


Bono et al. 

<: 

003 

) 

Catelan et al. 

(2004) 


Theoretical or semi-theoretical relations 


-2.101 0.231 +0.012 -0.770 + 0.044 

-2.353 0.175 -0.869 


-0.58 ± 0.04 


Dall’Ora et al. 

d200 

4) 

Del Principe et al. ( 

2006 

Sollima et al. ( 

2006 

i 3 


Sollima et al. ( 

2008 

i 5 


Borissova et al 

.(2009 

,3.7 


This paper (LMC) 1 7 

This paper (MW) 8 


-2.16 + 0.09 
-2.71+0.12 
-2.38 ± 0.04 
-2.38 ± 0.04 
-2.11+0.17 
-2.73 ± 0.25 
-2.53 ± 0.36 


Empirical relations 


0.12 + 0.04 
0.08 + 0.11 
0.08 + 0.11 
0.05 ± 0.07 
0.03 ± 0.07 
0.07 ± 0.04 


-1.05 + 0.13 
-1.07 + 0.11 
-1.05 

-1.06 + 0.01 
-0.95 + 0.14 


-0.56 ± 0.02 
-0.57 ± 0.02 
-0.57 ± 0.03 6 
-0.56 ± 0.02 
-0.56 ± 0.03 
-0.55 ± 0.06 9 
-0.56 ± 0.06 9 


'Zero-point of the original relation from the literature in the form: alogP + /;[Fe/H | + ZP 

2 Zero-point of the relation in the form: «(logP + 0.28) + Z?([Fe/H] + 1.58) + ZP, as recalibrated by Benedict et al. (2011 1 
3 Near-infrared photometry in the K s band. 

4 Metallicity is on the Zinn & West (1984) metallicity scale 


5 Metallicity is on the Carretta & Gratton (1997 metallicity scale 
6 Zero-point was calibrated by 


Benedict et al. 


(2011) neglecting the metallicity term 


7 Metallicity on the scale adopted by Gratton et al. (2004), which is , on average, 0.06 dex higher than Zinn & West (1984) scale 
8 Metallicities from < 


Clementini et al. 


(1995 i, they are on the scale of the high dispersion spectra (i.e., the 


Carretta et al. 


2009) 


9 Zero-point was calibrated by us considering only four RR Lyrae stars (XZ Cyg, UV Oct, SU Dra and RR Lyr) from 
|et al.|(|20TT) and excluding RZ Cep, since there are concerns about the parallax of this star 


Benedict 
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4.1. Method 


Using parallax data for calibration of a PL relation is complicated by the presence of statistical biases 
(e.g. Lutz & Kelker|1973] l and sample selection effects (e.g. |Malmquis t|1936 1. Non-linear transformations 
on the parallax cause a highly asymmetric uncertainty on the absolute magnitude when calculated using 
parallax and apparent magnitude information via the relation: m - M =-5 - 5 log (nr), where m is the 
parallax. Additionally, stars with a negative parallax measurement can not be used to calculate an absolute 
magnitude, though they do contain information. For these reasons, calculating an absolute magnitude for 
each star and fitting a PL relation directly to period and absolute magnitude leads to a biased result. 


An unbiased solution can be achieved through modelling the data and inferring the slope and zero-point 
of the relation via statistical methods. For a catalogue of N stars we can define x = (x\,X 2 , —,Xn) where 
the vector (jc, = m, 1, b, P, nr, A) describes the observed data on each object. P is the period, m the apparent 
magnitude, / and b the position, m the parallax, and A the extinction). We can additionally define that the 
vector (jco = mo, Iq, bo, P( h tq, Aq) gives the ‘true’ underlying object properties. 


We assume that the stars follow a PL relation of the form Mq = plogPo+d, although this can be changed 
to include other terms, such as metallicity, as needed. We can therefore model the true absolute magnitudes 
of the population as being Normally distributed around the PL relation, with the dispersion describing the 
intrinsic scatter on the relation: 

1 Q5 / M 0 -(plogj»o + « \ 2 

<Pm(x 0 \p, 6, crp L ) = -— e A ,TpL > (9) 

o~pi y2n 

where cr PL is the intrinsic dispersion of the PL relation. The parameters p and 6 are the slope and zero-point 
of the PL relation, which are to be found. 


The true absolute magnitude is calculated through: 


M 0 = m 0 + 5 log(tDo) + 5 - A 0 


( 10 ) 


The observations are Normally distributed around the true values with a standard deviation given by 
the formal error on the measurement: 


£(*1*0) = 


1 

e m e m e A ( 2 n) 3/2 



m—niQ 

em 



S(k,bo,Po) 


( 11 ) 


Assuming negligible errors on the position and period, the observations are described by a delta function. 
The terms e m , e m , ex are the formal errors on the parallax, magnitude and extinction. 


With the above models defined, the joint probability density function for the observations is: 

P(xi\p, 6, (tpl) - S(x) I <p M (xo\p,8,o-p L )S(x\xo)dxo (12) 

Jvx o 

the ‘true’ parameters xo are never known and so these values are marginalised through integration. The 
term S(x) is the selection function, which takes the probability of observing a star into account, given the 
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properties of the star and the instrument’s observational capabilities. To take the fact that Gaia is a magnitude 
limited sample into account, a step function is used with 


S(x) = 


0 , 


if G<20. 
otherwise. 


(13) 


The Maximum Likelihood Estimation of the parameters are found by maximising equation[l2]by vary¬ 
ing the parameters (a, p, ctpl)- This formulation avoids non-linear transformations on error effected data, 
and includes a selection function which avoids the Malmquist bias. 


4.2. Simulated Gaia data 


In order to check the application of the method defined in Section 4.1 we have used the sample of 23 
RR Lyrae stars in the MW discussed in Section |331 (see also Table[2]). In order to investigate the performance 
of the Gaia satellite and the contents of the end-of-mission catalogue, Gaia’s Data Processing and Analysis 
Consortium (DPAC) has a group working on the simulation of several aspects of the Gaia mission. One 


major product of this work is the Gaia Object Generator (GOG; Luri et al 2014), designed to simulate 
both individual Gaia observations and the full contents of the end-of-mission catalogue. GOG includes a 
full mathematical description of the nominal performance of the Gaia satellite, and is therefore capable of 
determining the expected precision in astrometric, photometric and spectroscopic observations. In general, 
the precision depends on the apparent magnitude of the star, its colour, and its sky position, which affects 
the number and type of observations made (due to the Gaia scanning law). 


To obtain a distance for each RR Lyrae star from the sample, we use: 


M k = -2.53 logP - 0.95 


(14) 


as determined in Equation[7]to obtain an absolute magnitude (neglecting the metallicity term for simplicity). 
We then determine a distance by combining this absolute magnitude with the apparent magnitude and ex¬ 
tinction as defined above. Colour information as (V-I) is obtained from the Hipparcos catalogue (Perryman 
and ESAj 1997[ ) where available. The apparent magnitude, position, colour, and period data form the basis 
of a synthetic catalogue of RR Lyrae stars, along with the distance obtained from the PMkZ relation, and is 
used as the input catalogue of ‘true’ parameters for GOG. 

GOG then creates simulated Gaia observations for our sample. We take the PMkZ elation (Equation[T4]) 
as true, as a study of the possible precision in PMkZ calibration after the Gaia data will become available. 


Using the fitting method described in Section 4.1 to the data including the simulated parallax observa¬ 
tions and simulated errors applied to parallax and apparent magnitude, we find a PMk relation of: 


M k = (-2.531 ± 0.038) log P + (-0.95 ± 0.01) 


(15) 
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Comparison of these results to the input relation shows very good agreement. It proves that the fitting 
procedure given in Section 4.1 is accurate and unbiased. When Gaia parallaxes will become available for 
the much larger sample of RR Lyrae stars, we will apply the described method to fit the PLkZ relation of 
RR Lyrae variables in the MW. Moreover, precise distance to the LMC obtained from the combination of 
Gaia parallaxes for a large sample of the bright LMC stars, will allow us to calibrate zero-point of the PLk s Z 
relation based on the LMC RR Lyrae stars. This provides a flavor of what will be possible to achieve with 
Gaia parallaxes. 


5. Summary 


We studied a sample of 70 RR Lyrae stars in the LMC, for which multi-epoch K s photometry from the 
VMC survey, precise periods from the OGLE III catalogue and spectroscopically determined metallicities 
(Gratton et a l.f|2004 >, are available. There are 13 epoch data in the K s band for all stars in the sample, that 
allowed us to determine mean K s magnitudes with a great accuracy. 


Specifically for this work we developed a fitting approach. This method has several advantages com¬ 
pared to the Minimum-Least-Squares fitting, such as taking into account potentially significant intrinsic 
dispersion of the data, non-negligible errors in two dimensions and the possibility of inaccuracies in the 
formal error estimates. We used this method to derive the PLk s Z relation of the 70 RR Lyrae stars in the 
LMC. Potentially the method could be used to fit any other sample of data. 


The zero-point of the derived PLk s Z relation was estimated in two different ways: (i) by assuming the 
distance to the LMC determined by Pietrzynski et al. (20131; (ii) by applying HST parallaxes of four MW 
RR Lyrae stars by Benedict et al. ( j201l| . The zero-point derived using the MW RR Lyrae stars is 0.2 mag 
larger and, consequently, gives a longer distance to the LMC: (m - M) = 18.68 ± 0.10 mag. In future studies 
we suggest to use the relation based on the precise distance to the LMC: 


M Ks = (-2.73 ± 0.25)logP + (0.03 ± 0.07)[Fe/H] 

- (1.06 ±0.01) (16) 

We found a negligible dependence of the Mk s on metallicity, which could be caused by the relatively 
small range in metallicity covered by the LMC RR Lyrae stars. Thus, we applied the fitting approach to 23 
RR Lyrae stars in the MW, for which absolute Mk and My magnitudes are known from Baade-Wesselink 
studies. We derived the PLk s Z relation for MW RR Lyrae stars in the form: 


M k = 


(-2.53 ± 0.36)logP + (0.07 ± 0.04)[Fe/H] 
(0.95 ±0.14) 


(17) 
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Even though the metallicities of the MW RR Lyrae stars span a wide range [-2.5; 0.17] dex, the de¬ 
pendence on metallicity is relatively small and consistent, within the errors, with the slope in metallicity 
found for the LMC RR Lyrae variables. We concluded that the small range of metallicities doesn’t cause the 
negligible dependence of the Mk on metallicity for the LMC RR Lyrae stars. 

To solve the problem of the PLk s Z zero-point, a large sample of RR Lyrae stars with precisely de¬ 
termined parallaxes is necessary. A great contribution to this field is expected by the Gaia satellite. By 
using GOG we simulated Gaia parallaxes of 23 MW RR Lyrae stars with observational errors. We present a 
method for the calibration of the PL relation which avoids several of the problems which arise when using 
parallax data. The method was tested by deriving the PLk s relation based on the simulated Gaia parallaxes. 
When combined with metallicity and photometry from other sources and a statistical tool such as the one 
developed in the present study, the extraordinary large sample of Gaia parallaxes for RR Lyrae stars will 
allow us to estimate these relations with unprecedented precision. 
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A. Bayesian fitting approach 


This method is based on the prescription of Hogg et al. (2010), extended into three dimensions and 
implemented in Python using a Markov Chain Monte Carlo sampler to obtain parameter estimates along 
with their complete posterior PDF. Initially, we model the data as being drawn from a thin plane defined by: 


f(x, y)=Ax + By + C 


(Al) 


where A is the slope in the x axis, B is the slope in the y axis, and C is the intercept. In this initial model we 
assume that we have data in three axis, x, y and z, with errors only in the z axis. 


In this model, given an independent position ( Xi,y : - ), an uncertainty <r Z i, slopes A and B, and an intercept 
C, the frequency distribution p(zi\xi, //,, <r z i,A, B, C) for Zi is 

, , „ 1 [Zi-Axi-Byi-C] 2 ) 

p(Zi\x{, y h cr zi ,A, B, C) = exp - — 5 - (A2) 


yl 2ncr2 zi 


2 a 1 . 

Zl 
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Therefore the likelihood is defined as: 


N 

1=1 


(A3) 


Taking the logarithm, 


N 

In if = K-^j 

i=i 


[Zi - A Xj-B iji - C]- 
2 cr 2 . 

ZI 


(A4) 


which is effectively the least-squares solution. K is a normalisation coefficient. Returning to Bayes rule it is 
possible to define: 

p({z i }l l \A,B,C,I)p(A,B,C\I) 


p(A, B, C\{zi}f =l , T) = 


p({Zi}? =l \I) 


(A5) 


ki)f = i is all the data Zi- / is all of the information of x and y, (x,-, yi\^ =] , along with the formal errors 
[cr Xi , ( T yi , ovl'^p plus any other prior information which may be available. In our case we use uninformative 
(uniform) priors, making our inference method analogous to Maximum Likelihood Estimation. 


A.l. Multiple errors, no dispersion 


As in this case there exist errors in more than one axis, they can be put together into a covariance tensor 


Si 


9 



cr. 

XI 

0~xyi 

o-xzi 

0~xyi 

cr 2 . 

yt 

°~yzi 

°~xzi 

o'yzi 

cr 2 

zi 


(A6) 


With errors in several dimensions, our observed data point ( Xi,yuZi) could have been drawn from any true 
point along the plane ( x,y,z ). Making the probability of the data, given the model and the true position: 


p{xi,iji,Zi\Si,x,y,z) 


1 

- , exp 

2 :r VdettS,) 


-[ Zi-ZfSf 1 |Z,-Z| 


where we have implicitly made column vectors 


(A7) 



X 


Xi 

Z = 

y 

; Z; = 

Ui 


z 


. Zi 


(AS) 


In only two dimensions (e.g. y and z), the slope (e.g. B) can be described by a unit vector v orthogonal to 
the line or linear relation (at any x): 


1 

V = - 

-B 


- sind 

Vl +B 2 

1 


cos 6 


(A9) 
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where the angle 9 = arctan B is made between the line and the y axis. The orthogonal displacement A,- of 
each data point (j/,-, zi) from the line is given by 


A; 


= v 


y> 

Zi 


- C cos 6 


(A10) 


Instead of extending fully into three dimensions, we will assume a negligible error in x (which will be the 
period, so has justifiably higher precision). The value of x can therefore be input directly into A; without 
worrying about the interplay between the other parameters. 


A; = V 1 


Vi 

Zi 


- (C cos 6 + A xi) 


(All) 


Assuming negligable errors in x also redefines the covariance matrix of the errors as: 


Si = 


yi °~yu 

O-y-J (T 2 zi 


(A 12) 


Similarly, each data point’s covariance matrix S, projects down to an orthogonal variance I? given by 


if = v T Sii> 


(A 13) 


and then the log likelihood for {A, B, C) or ( A,9,C cos 9) can be written as 


lnJz? - K - 


N 

z 

1=1 


21? 


where K is some constant. This likelihood can be maximized to find A, B and C. 


(A 14) 


A.2. Dispersion 


The final step is to introduce an intrinsic variance in the line, V, orthogonal to the line. 


According to Hogg et al. (2010), each data point can be treated as being drawn from a projected 
distribution function that is a convolution of the projected uncertainty Gaussian, of variance I? defined 
above, with the Gaussian intrinsic scatter of variance V. Therefore the likelihood becomes: 


N i N 

InJSf = K - ^ - ln(S? + V) - £ 


A 2 


1=1 


i= 1 


2 [I? + V] 


(A 15) 


where again K is a constant, everything else is defined as above. We then solve for A, B , C, and V by 
maximising the log likelihood. The optimisation is performed using the adaptive MCMC sampler EMCEE 
developed by Foreman-Mackey et al. ( 2013| ). Any optimisation algorithm (e.g. Nelder-Mead, Powell, etc.) 
will find the maximum of the log likelihood. MCMC was chosen due to the evaluation of the full posterior 
PDF of the parameters, which is useful for the determination of formal errors. 
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B. Metallicities for the MW RR Lyrae stars analyzed with the B-W method 

In this Appendix we provide a summary of spectroscopic metal abundances ([Fe/H]) derived for the 
field RR Lyrae stars analyzed with the B-W technique, by studies mainly based on high resolution spectro¬ 
scopic material. Exceptions are the values with reference 1, 2 and 4 that come from compilations which 
include metallicities measured from low resolution spectroscopic data and photometric indices (see discus¬ 
sion in Section 4.2 of Cacciari et al. 1992). 

By applying our fitting approach we derived the PMkZ relation for 23 MW RR Lyrae stars described 
in Section [33] adopting as an alternative the metallicity values from Lambert et al. (19961: 


Lambert et al. (1996 


M k = (-2.66 ± 0.36)logP + (0.05 ± 0.04)[Fe/H] 
- (1.00 ±0.15) 


(Bl) 


The intrinsic dispersion of the relation is found to be 0.007 mag. The RMS deviation of the data around 
the relation, neglecting the intrinsic dispersion, is 0.090 mag. The slope in logP of the relation based on 
Clementini et aL|( jl995 i’s metallicities (Eq. [7]) differs from the slope obtained using Lambert et al. ( 1996j )’s 


metallicities (Eq. B1), however the values are consistent within the errors. 
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Table 4. Literature metallicities of the MW RR Lyrae stars 
analyzed with the B-W method 


Star 

[Fe/H] 

O' [Fe/H] 

Reference 


(dex) 

(dex) 



UU Cet 


SW And 


RR Cet 


X Ari 


-1.38 

0.08 

(5), average of Fel and Fell 

-1.33 

0.08 

(3), average of Fel and Fell 

-0.95 

0.16 

(7*), re-calibration of AS and As value adopted by (7) 

-1.45 

0.16 

(7*), re-calibration of AS in (7) and As value adopted by (5) 

-1.20 

0.20 

(4) 

-1.28 


(1.2) 

-0.06 

0.08 

(5), average of Fel and Fell 

-0.41 

0.10 

from Fel in (7) 

-0.24 

0.12 

from Fell in (7) 

-0.34 

0.085 (s=0.12) 

average of Fel and Fell in (7) 

-0.27 

0.15 

(11), from Fel and solar abundance 7.51 

-0.24 

0.15 

(11), from Fell and solar abundance 7.51 

-0.255 

0.15 

(11). average of Fel and Fell 

+0.20 

0.08 

(3), average of Fel and Fell 

-0.16 


(9) 

-0.15 

0.15 

(4) 

-0.24 


(1.2) 

-1.38 

0.08 

(5), average of Fel and Fell 

-1.62 

0.15 

(11), from Fel 

-1.49 

0.15 

(11), from Fell 

-1.18 

0.09 

(3), average of Fel and Fell 

-1.61 


(9) 

-1.36 

0.16 

(7*), re-calibration of AS and As value adopted by (7) 

-1.25 

0.10 

(4) 

-1.45 


(1.2) 

-2.50 

0.08 

(5), average of Fel and Fell 

-2.19 

0.17 

(6), from Fel 

-2.54 

0.09 

from Fel in (7) 

-2.75 

0.08 

from Fell in (7) 

-2.66 

0.105 (s=0.15) 

average of Fel and Fell in (7) 

-2.74 

0.09 

(3), average of Fel and Fell 

-2.61 


(8), from Fel 

-2.62 


(8), from Fell 

-2.68 


(9) 

-2.20 

0.10 

(4) 

-2.43 


(1.2) 



Table 4—Continued 


Star 

[Fe/H] 

<T [Fe/H] 

Reference 


(dex) 

(dex) 



AR Per 

-0.34 

0.16 

(5*), re-calibration of AS and As value adopted in (5) 


-0.23 

0.07 

from Fel in (7) 


-0.41 

0.08 

from Fell in (7) 


-0.31 

0.09 (s=0.13) 

average of Fel and Fell in (7) 


-0.24 

0.15 

(11) from Fel 


-0.29 

0.15 

(11) from Fell 


-0.32 


(9) 


-0.30 

0.20 

(4) 


-0.30 


(L2) 

RX Eli 

-1.63 

0.16 

(5*), re-calibration of AS and AS value adopted in (5) 


-1.98 

0.16 

(7*), re-calibration of AS and AS value adopted in (7) 


-1.40 

0.20 

(4) 


-1.33 


(L2) 

RR Gem 

-0.32 

0.16 

(5*), re-calibration of AS and AS value adopted in (5) 


-0.44 

0.16 

(7*), re-calibration of AS and AS value adopted in (7) 


-0.30 

0.25 

(4) 


-0.29 


(1,2) 

TTLyn 

-1.64 

0.16 

(5*), re-calibration of AS and AS value adopted in (5) 


-1.63 

0.08 

from Fel in (7) 


-1.33 

0.06 

from Fell in (7) 


-1.44 

0.150 (s=0.21) 

average of Fel and Fell in (7) 


-1.64 

0.15 

(11), from Fel 


-1.53 

0.15 

(11), from Fell 


-1.41 


(9) 


-1.35 

0.20 

(4) 


-1.56 


(1,2) 

TSex 

-1.20 

0.16 

(5*), re-calibration of AS and AS in (5) 


-1.75 

0.12 

from Fel in (7) 


-1.50 

0.09 

from Fell in (7) 


-1.59 

0.125 (s=0.18) 

average of Fel and Fell in (7) 


-1.20 

0.15 

(4) 


-1.34 


(1,2) 

RR Leo 

-1.37 

0.16 

(5*), re-calibration of AS and AS value in (5) 
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Table 4—Continued 


Star 

[Fe/H] 

O' [Fe/H] 

Reference 


(dex) 

(dex) 



WY Ant 


WCrt 


TU UMa 


UU Vir 


-1.54 

0.11 

from Fel in (7) 

-1.17 

0.10 

from Fell in (7) 

-1.34 

0.185 (s=0.26) 

average of Fel and Fell in (7) 

-1.39 


(9) 

-1.15 

0.20 

(4) 

-1.60 


(1,2) 

-1.32 

0.16 

(5*), re-calibration of AS and AS value in (5) 

-1.96 

0.10 

(10), from Fel 

-1.96 

0.10 

(10), from Fell 

-1.96 

0.10 

average of Fel and Fell in (10) + errors front us 

-1.36 

0.16 

(7*), re-calibration of AS and AS value in (7) 

-1.25 

0.20 

(4*) 

-1.48 


(1,2) 

-0.89 

0.16 

(5*), re-calibration of AS and AS value in (5) 

-0.91 

0.16 

(7*), re-calibration of AS and AS value in (7) 

-0.70 

0.20 

(4*) 

-0.54 


(1,2) 

-1.38 

0.16 

(5*), re-calibration of AS and AS value in (5) 

-1.64 

0.08 

from Fel in (7) 

-1.45 

0.08 

from Fell in (7) 

-1.55 

0.095 (s=0.13) 

average of Fel and Fell in (7) 

-1.31 

0.05 or 0.14 

(6), from Fel, error of 0.05 likely a typo 

-1.72 

0.15 

(11), from Fel 

-1.57 

0.15 

(11), from Fell 

-1.46 


(9) 

-1.25 

0.20 

(4) 

-1.51 


(1,2) 

-0.64 

0.16 

(5*), re-calibration of AS and AS value in (5) 

-0.85 

0.12 

from Fel in (7) 

-0.79 

0.07 

from Fell in (7) 

-0.81 

0.03 (s=0.04) 

average of Fel and Fell in (7) 

-0.90 


(9) 

-0.50 

0.15 

(4) 

-0.87 


(1,2) 
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Table 4—Continued 


Star 

[Fe/H] 

(dex) 

C [Fe/H] 

(dex) 

Reference 

SW Dra 

-0.91 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-1.37 

0.15 

(11), from Fel 


-1.28 

0.15 

(11), from Fell 


-0.81 

0.16 

(7*), re-calibration of AS and AS value in (7) 


-1.14 

0.20 

(4) 


-1.12 


(1,2) 

RV Oct 

-1.92 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-1.54 

0.10 

from Fel in (10) errors from us 


-1.54 

0.11 

from Fell in (10) errors from us 


-1.54 

0.11 

average of Fel and Fell in (10) errors from us 


-1.98 

0.16 

(7*), re-calibration of AS in (7) and AS value 


-1.75 

0.20 

(4*) 


-1.71 


(1,2) 

TV Boo 

-2.31 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-2.36 

0.16 

(7*), re-calibration of AS in (7) and AS value 


-2.30 

0.15 

(4) 


-2.44 


(1,2) 

RS Boo 

-0.37 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-0.55 

0.11 

from Fel in (7) 


-0.39 

0.08 

from Fell in (7) 


-0.45 

0.08 (s=0.11) 

average of Fel and Fell in (7) 


-0.48 

0.15 

(11), from Fel 


-0.33 

0.15 

(11), from Fell 


-0.40 

0.25 

(4) 


-0.36 


(1,2) 

VY Ser 

-1.71 

0.08 

(5), average of Fel and Fell 


-2.09 

0.08 

from Fel in (7) 


-1.76 

0.06 

from Fell in (7) 


-1.88 

0.165 (s=0.23) 

average of Fel and Fell in (7) 


-1.71 

0.07 

(3), average of Fel and Fell 


-2.00 

0.15 

(11), from Fel 


-1.90 

0.15 

(11), from Fell 


-1.80 

0.15 

(4) 


-1.79 


(1,2) 
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Table 4—Continued 


Star 

[Fe/H] 

(dex) 

O" [Fe/HJ 

(dex) 

Reference 

V445 Oph 

+0.17 

0.08 

(5), average of Fel and Fell 


+0.13 

0.10 

(3), average of Fel and Fell 


-0.26 

0.16 

(7*), re-calibration of AS in (7) and AS value in (7) 


+0.24 


(9) 


-0.30 

0.25 

(4) 


-0.19 


( 1 . 2 ) 

TW Her 

-0.58 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-0.56 

0.16 

(7*), re-calibration of AS and AS value in (7) 


-0.50 

0.15 

(4) 


-0.69 


( 1 , 2 ) 

AV Peg 

-0.03 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-0.35 

0.09 

from Fel in (7) 


-0.04 

0.06 

from Fell in (7) 


-0.14 

0.155 (s=0.22) 

average of Fel and Fell in (7) 


0.00 

0.20 

(4) 


-0.08 


( 1 . 2 ) 

RVPhe 

-1.69 

0.16 

(5*), re-calibration of AS and AS value in (5) 


-1.75 

0.16 

(7*), re-calibration of AS and AS value in (7) 


-1.35 

0.25 

(4) 


-1.69 


( 1 . 2 ) 


Note. — (1) [Fe/H] from Table 1 in |Femley et al.|(T998a} ; (2) [Fe/H] values from 
|Feast et aT7| ([2008 ); (3) [Fe/H] values from |Nemec et al.|(20 1 3| using the values from 
the VWA analysis as they are listed in column (9) of Table 7 in that paper. They are 
the average of the Fel and Fell abundances. (4) [Fe/H] values from Table 11 of 
ICacciari et al~| fj~ 992 > ; ( 4 *) [p e /H] values from Table 16 of |Skillen et ah] ( fl993l >; 
(5) [Fe/H] values from abundance analysis of high resolution spectra performed by 
|Clementini et a 1. | ( | 1995[ >. Values are the average of the Fel and Fell measurements 
adopting for the solar abundance: log e(FeI)=7.56 and log e(FeII)=7.50, respec¬ 
tively (see Table 12 of |Clementini et al.|1995j ); (5*) [Fe/H] values obtained from 
|Clementini et a 1. | ( | 1995[ re-calibration of the AS index. ( 6 ) [Fe/H] values from Ta¬ 
ble 7 of [Pancino et al.||2015| l. Abundances are from Fel averaging values from 
different spectra of the same star as described at the end of the paper Section 4.4. 
The error for the abundance of TU Uma is likely a typo, according to the paper Ta¬ 
ble 6 likely should be 0.14 dex; (7) Metal abundances from|Lambert et aT7| (j 1996[. 
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[Fe/H] values were derived from the photometric determinations in the paper Ta¬ 
ble 3 adopting for the sun log e(Fe)= 7.51 (according to what stated in the footnotes 
of the paper Table 5); (7*) [Fe/H] values obtained from |Lambert et ah] l | 1996) re¬ 
calibration of the AS index (Equation 3 in that paper) which was derived by these 
Authors using the Fell abundances and AS from |Blanco| ( j 1992| . [Lambert et al.| 
( | 1996| > does not provide errors for the metallicities from AS, hence we adopted an 
error of 0.16 dex, as done by Clementini et ak|l |1995[ l for their metallicities from AS; 
(8) [Fe/H] values from Table 10 of |Haschke et al. | ( [20 1 2[ ; (9) [Fe/H] values abun¬ 
dance analysis performed by Wallerstein & Huang (2010), no errors are provided; 
(10) [Fe/H] values from |For et al.|{20TT| obtained as the average weighted by errors 
of the values in the paper Table 5. Errors are the sum in quadrature of the individual 
errors in Table 5 divided by the square root of N (with N number of measurements, 
i.e.: 11 for WY Ant and 17 for RV Oct). The average value from Fel for WY Ant 
published in Table 11 of |For et aT]j20 11 is -1.95 dex while we find -1.96 dex with 
our procedure. We also list the abundances from Fel and Fell separately, with errors 
calculated as per the above procedure; (11) [Fe/H] values from |Femley & Bames| 
( | 1996| >.They are from Fell with values taken from the paper Table 4b. Those from 
Fel are taken from the paper Table 4a. Errors are estimated by the authors to be of 
±0.15 dex in [Fe/H]. 
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